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Abstract: With realizing the importance of ecosystem services to society, the efforts to evaluate the 
ecosystem services have increased. As the largest tributary of the Yellow River, the Weihe River has been 
endowed with many ecological service functions. Among which, water yield can be a measure of local 
availability of water and an index for evaluating the conservation function of the region. This study aimed 
to explore the temporal and spatial variation of water yield and its influencing factors in the Weihe River 
Basin (WRB), and provide basis for formulating reasonable water resources utilization schemes. Based on 
the InVEST (integrated valuation of ecosystem services and tradeoffs) model, this study simulated the 
water yield in the WRB from 1985 to 2019, and discussed the impacts of climatic factors and land use 
change on water yield by spatial autocorrelation analysis and scenario analysis methods. The results 
showed that there was a slight increasing trend in water yield in the WRB over the study period with the 
increasing rate of 4.84 mm/10a and an average depth of 83.14 mm. The main water-producing areas were 
concentrated along the mainstream of the Weihe River and in the southern basin. Changes in water yield 


were comprehensively affected by climate and underlying surface factors. Precipitation was the main factor 
affecting water yield, which was consistent with water yield in time. And there existed significant spatial 
agglomeration between water yield and precipitation. Land use had little impact on the amount of water 
yield, but had an impact on its spatial distribution. Water yield was higher in areas with wide distribution 
of construction land and grassland. Water yield of different land use types were different. Unused land 
showed the largest water yield capacity, whereas grassland and farmland contributed most to the total 
water yield. The increasing water yield in the basin indicates an enhanced water supply service function of 
the ecosystem. These results are of great significance to the water resources management of the WRB. 
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1 Introduction 


Ecosystems directly or indirectly provide various services to human beings (Symmank et al., 
2020), which were defined as "the benefits that humans derive from nature" (Hassan and Scholes, 
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2005). As the attribute of the ecosystem and the key component of an ecological process, water 
yield can produce water-related ecosystem services (Costanza et al., 1998). For example, it 
regulates surface runoff over the dry and flood seasons, reduces potential drought and flood risks, 
and ensures drinking water sources (Zhang et al., 2010). In summary, water yield is a resource 
supporting human life and development, such as agriculture, industry, human consumption, 
hydropower, fisheries, and recreational activities (Natalia et al., 2015). 

The growth in the global human population, the improvement in living standards, changes in 
consumption patterns, and the expansion of irrigated agriculture over the last few decades has 
resulted in a gradual increase in the demand for water resources (Ercin and Hoekstra, 2014; 
Gomez et al., 2014; Vorosmarty et al., 2000). Consequently, many countries are facing increasing 
water scarcity challenges (Qiu, 2010; Liu and Wu, 2012; Mekonnen and Hoekstra, 2016). In 
China, high demand for food and the imbalance of regional economic development have resulted 
in increased water resources utilization and uneven distribution (Mekonnen and Hoekstra, 2016). 
Relevant studies have shown a sharp decrease in measured runoff in major river basins of China 
in recent decades (Xu et al., 2010). The changes to water resources of the Yellow River have 
attracted more attention for the importance of this river in China (Qiang et al., 2009). In the late 
1990s, the average runoff entering the Yellow River was only 35.99x10° m°, a decrease of 
62.06x10° m?, compared with the 1950s (Xia et al., 2007). And this runoff further decreased in 
the 21* century (Zhao et al., 2018). 

The Weihe River is the largest tributary of the Yellow River and provides irrigation water for 
nearly 1x10* km? of farmland in Guanzhong Plain, supporting nearly 61% of the population, 56% 
of cultivated land, and 81% of gross domestic product (GDP) in Shaanxi Province, China (Liu 
and Hu, 2008). However, the imbalance between water supply and demand in the Weihe River 
Basin (WRB) have intensified due to impacts of rapid human development, the construction of 
water conservancy, soil conservation projects, and other engineering measures (Ma et al., 2008; 
Cheng et al., 2019; Xu et al., 2021). Information in the Shannxi Water Resources Bulletin (2019) 
states that agriculture is the sector consuming the most water in the basin (55.13x10° m°), which 
includes forestry, animal husbandry, fishery, and livestock farming, accounting for 60% of total 
water consumption; The water consumption of industry and municipal use was 14.85x108 and 
14.28x108 m3, accounting for 16.05% and 14.12% of total water consumption, respectively. Water 
consumption in the basin will increase further with the development of economy, imposing more 
pressure on available water resources (Zhang et al., 2016). 

The WRB belongs to a typical transition climate zone in the arid and semi-arid area of 
northwestern China, and water resources in this area are sensitive to climate change (Chen et al., 
2013). Global climate variability over the last few centuries has been characterized by 
temperature rise and precipitation changes. These changes have had a destructive impact on the 
natural ecosystem, human economic, and social development (Milliman et al., 2008). Water yield 
has gradually become an important factor restricting the sustainable development of society. 
Therefore, an evaluation of water yield of the WRB is of great significance. 

Model can be a useful tool for the evaluation of water-related ecosystem service functions and 
the quantitative estimation of water yield under various conditions. The models include soil and 
water assessment tool (SWAT) (Baker and Miller, 2013; Gassman et al., 2017), social values for 
ecosystem services (SolVES) (Greg et al., 2012; Sherrouse et al., 2011), multi-scale integrated 
models of ecosystem services (MIMES) (Boumans and Costanza, 2008), and integrated valuation 
of ecosystem services and tradeoffs (InVEST) (Sharp et al., 2020). Among these models, the 
InVEST model has been the most widely applied due to its less data requirement and ability to 
visualize simulation results (Huang et al., 2013; Scordo et al., 2018; Cara et al., 2020). The 
InVEST model integrates a series of ecosystem processes, and can simulate the quality and value 
of ecosystem services regulated by land use, physical environmental factors, and socio-economic 
factors (Jiang et al., 2021). Numerous studies have applied it to various regions of the world and 
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have achieved good results, such as North Korea (Kim and Jung, 2020), the Wildcat Creek 
Watershed in Indiana, the Upper Upatoi Creek Watershed in Georgia, the United States 
(Dennedy-Frank et al., 2016), and some regions of China (Yang et al., 2019; Yin et al., 2020; Li et 
al., 2021). However, there have been few attempts to apply the InVEST model to the WRB. Those 
few studies that have applied the InVEST model to the WRB were not able to identify the 
temporal and spatial variation characteristics in water yield due to their short study periods. Most 
of these studies also only analyzed the quantitative relationships between water yield and driving 
factors (Yang et al., 2019), while not analyzing their spatial correlations. 

The present study used the InVEST model to (1) simulate water yield of the WRB from 1980 
to 2019, and analyze its spatiotemporal variation; (2) identify the main factors regulating water 
yield under different climate factors and land use types; and (3) discuss the mechanism of 
regulating water yield. 


2 Materials and methods 
2.1 Study area 


The Weihe River originates from Niaoshu Mountain, Weiyuan County, Gansu Province, following 
which it flows through the Ningxia Hui Autonomous Region and Shaanxi Province from east to 
west, including 10 regions and 84 counties. The Weihe River merges into the Yellow River 
eventually at Tongguan County, Shaanxi Province, China. As the largest tributary of the Yellow 
River, the Weihe River has a total length of 818 km, and drains an area of 13.5x10* km?. WRB 
located in the southeastern Loess Plateau, China (33°42'—37°20'N, 106°18’-110°37’E). The 
terrain of the WRB gradually decreases from west to east. And there is an elevation difference of 
over 3000 m a.s.l. (Fig. 1). The climate of the WRB is cold and dry in winter, and hot and rainy in 
summer. The annual average temperature and average annual precipitation are 7.8°C—13.5°C and 
572 mm, respectively. 
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Fig.1 Elevation, and hydrological and meteorological stations in the Weihe River Basin, China 
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2.2 Data sources 


The data required by the model include meteorological, soil and land use/land cover (LULC) data. 
The raster input data for the model were derived directly or indirectly from these data. Table 1 
summarizes the relevant basic data used in the present study. 


Table 1 Relevant basic data sources and description 


Data 


Data description 


Data source 


Climate data 


Soil data 


Land use/land cover 


Digital elevation model 
(DEM) 


Streamflow data 


Runoff data 


Daily precipitation 

Daily maximum temperature 
Daily minimum temperature 
Daily mean temperature 


Soil texture (clay, sand, and silt) 
Soil organic carbon 
Soil depth 


Land use/land cover during 1980-2020 at 1 km 
spatial resolution 


Elevation of the Weihe River Basin at 30 m 
spatial resolution 


Measured runoff of Huaxian and Zhuangtou 
hydrological stations from 1980 to 2019 


Restored runoff of Huaxian and Zhuangtou 


China Meteorological Science Data 
Sharing Service Network 
(http://www.cma. gov.cn/) 


Harmonized World Soil Database 
(HWSD) 


Resource and Environmental Science 
Data Center, Chinese Academy of 
Sciences (http://www.resdc.cn) 


Geospatial Data Cloud Platform of 
Chinese Academy of Sciences 
(https://www.gscloud.cn/) 


Yellow River Hydrological Yearbook 


Runoff data of the Yellow River Basin 


hydrological stations from 1980 to 1984 


2.3  InVEST model 


2.3.1 Introduction of the nVEST model 


Budyko curve (Budyko, 1974), and annual average precipitation (P) were used in the InVEST 
model to calculate water yield. Water yield simulated by the model represents the water outflow 
from the landscape, including surface flow, underground flow, and base flow (Sharp et al., 2020). 
The model assumes that water yield of a pixel reaches the specified outlet through one of several 
paths mentioned above (Zhang et al., 2004; Donohue et al., 2012). The formula is expressed as: 


AET 

Y(x)= ja SEUI x P(x), (1) 
P(x) 

where AET(x) is the annual actual evapotranspiration (mm); and P(x) is the precipitation for a 
pixel (mm). 


For vegetation-covered land, SO was calculated (Fu, 1981; Zhang et al., 2004): 
x 
a 1/@ 
AET(x) 4 PET(x) lee PET(x) l (2) 
P(x) P(x) P(x) 


where PET(x) is the potential evapotranspiration (mm); and w(x) is a nonphysical parameter 
characterizing regional climate and soil conditions, ranging from 1.25 to 5.00 (Yang et al., 2008; 
Donohue et al., 2012). The minimum value can be selected when the root depth is 0 cm (bare soil), 
and is calculated as: 


AWC(x) 
P(x) 
where Z is the local P and other hydrogeological features, with possible values of 1 to 30; and 

AWC(x) is the plant available water capacity (mm), and is calculated as: 


AWC(x) = min(rest layer depth, root depth) x PAWC , (4) 
where PAWC is the plant available water capacity (Fig. 2a); and rest layer depth is the depth of 


o(x)=Z +1.25, (3) 
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root restricted layer, often expressed as the depth at which 95% of root biomass of the plant (Fig. 
2b). 
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Fig. 2 Spatial distributions of biophysical characteristics of the Weihe River Basin, China. (a), PAWC, plant 
available water capacity; (b), soil depth; (c), land use/land cover (LULC) in 2020. 
2.3.2 Data preparation 
The data input to the model included rasterized annual P, average annual ETo, LULC, depth of the 
root restricting layer, PAWC, the layers of watersheds and sub-watersheds, the biophysical table, 
and an appropriate Z parameter (Sharp et al., 2020). All raster data used in the present study were 
unified into the Krasovsky_1940_Albers coordinate system before input. The sections below 
describe the processing methods used. 

(1) P and ETo 

The annual P and ETo data were obtained by adding the daily precipitation and monthly ETo 
data within a year, and Kriging was then used to interpolate the data of 22 meteorological stations 
around the WRB, which provided a value for each cell in the raster. ETo was calculated by the 
modified Hargreaves equation (Droogers and Allen, 2002): 


ET, =0.0013x 0.408 x RA x (T,, +17) x (TD —0.0123P)°”*, (5) 


where Tav is the average of mean Tmax and Tmin for every month (°C), respectively; TD is the 
difference between mean Tmax and Tmin for every month (°C), and RA is the extraterrestrial 
radiation (MJ/(m*d)). Radiation data were obtained from the United Nations Food and 
Agriculture Organization (FAO) Irrigation Drainage Paper (Fao et al., 1982). 

(2) PAWC 

PAWC (Fig. 2a) is the field capacity minus wilting point, ranging from 0 to 1, which was 
calculated as (Zhou et al.; 2005): 


PAWC = 54.509 — 0.132 x sand — 0.003 x (sand)? — 0.055 x silt — 0.006 x (silt)? — 


0.738 x clay + 0.007 x (clay)? — 2.688 x OM + 0.501x (OM)’, 


where sand, silt, clay, and organic matter (OM) were expressed as percentages (%). 

(3) LULC 

ArcGIS 10.2 was used to reclassify 25 secondary land use types into 6 primary land types (Fig. 
2c): farmland, forestland, grassland, water body, construction land, and unused land. 

(4) Watershed and sub-watershed 

In Arcgis10.2, through the hydrology tool, the WRB was finally divided into 377 
sub-watersheds by setting the threshold flow accumulation (TFA) to 4x105 m°. 

(5) Biophysical table 

The biophysical table (Table 2) contains the information for each LULC grid required by the 
model (Sharp et al., 2020). The LULC_veg determines which AET calculation formula to be used. 
LULC with vegetation cover is assigned as 1, and others are assigned as 0. K. is the 
evapotranspiration coefficient of vegetation, which is based on alfalfa, and used to adjust the ETo 
to obtain PET, with a range from 0.0 to 1.5. Root depth (Fig. 2b) is the maximum depth that plant 


(6) 
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root system can extend, which is often expressed as the depth at which 95% of root biomass 
occurs (Allen et al., 2006). The data in Table 2 were obtained from the relevant literature (Yang et 
al., 2020), and from the values recommended by the FAO and the InVEST model. 


Table 2 Biophysical table for the InVEST model 


Code Land use type LULC_veg K.: Root depth (m) 
1 Farmland 1 0.65 400 
2 Forestland 1 1.00 1500 
3 Grassland 1 0.65 500 
4 Water body 0 1.00 1 
5 Construction land 0 0.30 1 
6 Unused land 0 0.50 1 


Note: LULC_veg, land use/land cover with vegetation cover; Ke, evapotranspiration coefficient. 


2.3.3 Model calibration 


The InVEST model simulates natural runoff of a basin. The simulation accuracy largely depends 
on parameter Z (Eq. 3). The present study validated the parameter Z by comparing the simulated 
water yield data with natural runoff data at the Huaxian and Zhuangtou hydrological stations in 
the WRB from 1980 to 1984. The accuracy of simulation was assessed using the correlation 
coefficient (R?), Nash coefficient (NSE), and relative error (Re) (Rientjes et al., 2011). The study 
set the range of Z value for the study area to 3.6—9.0, consistent with those of relevant studies (Wu 
et al., 2018; Yue et al., 2021). The model was initiated with a Z value of 3.6, following which Z 
was gradually increased by 0.2. The process was continued until the maximized values of R?, NSE, 
and Re were obtained. After nearly a thousand tests, the optimal values of Re, R?, and NSE of 2%, 
0.98, and 0.63, respectively were obtained at an input Z value of 8.8. 


2.4 Trend analysis 


Mann-Kendall (M-K) test is a non-parametric method, which does not need samples to follow a 
certain distribution, and is not disturbed by a few outliers (Mann, 1945; Kendall, 1975). The M-K 
test is often used to identify trends in precipitation and drought under the influence of climate 
change (Sheng and Paul, 2004). This study applied M-K test to detect the trends in the water yield, 
and climate factors over time. It should be noted that the statistic Z value from the M-K test is 
different to the seasonal parameter Z required as input into the InVEST model. 


2.5 Correlation analysis 


2.5.1 Global spatial autocorrelation 

The present study used Moran's J to express global spatial autocorrelation. This index was used to 
analyze the overall correlations between spatial units, and to assess whether a spatial 
agglomeration existed. The valueof Moran's / ranges from —1 to 1. />0 indicates that spatial 
attribute value has a spatial agglomeration effect with the surrounding attribute value, and the 
closer to 1 the J value, the more significant the spatial agglomeration effect. /<O indicates that the 
spatial attribute value has a spatial difference effect with the surrounding attribute value. The 
spatial difference expressed by Moran's J becomes more significant, and the closer the value is to 
—1, whereas /=0 indicates no spatial autocorrelation (Tu and Xia, 2008). The spatial weight matrix 
was set as a simple binary adjacency matrix. The calculation of Moran's 7 is as follows (Moran, 


1950): 
jo lt DED IW, (x; -7)(x; =z) 
DiZi;  Di-z7 
where W; is the spatial weight matrix of cell i and j; x; and x; are the measured values of cell x and 
y, respectively; x is the average measured value of all cells; and n is the number of all 


, (7) 
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evaluation units. 

2.5.2 Bivariate local Moran's I 

We used local indicators of spatial association (LISA) method to assess the strength of the 
correlation between attributes of each spatial unit and the attributes of adjacent units according to 
local Moran's J (Harries, 2006). Here, 7;>0 indicates that spatial units are highly correlated, 
including high-high and low-low types; /;<O indicates large differences in spatial units, including 
high-low and low-high types. The formula is as follows a a et al., 2019): 


n? (%,-X) >) Wa- 


LOM; LIe- 

This method can be used to study the spatial agglomeration effect of meteorological factors and 

water yield. The high-high and low-low types refer to proportional relationships between 

meteorological factors and water yield, whereas high-low and low-high types indicate inverse 
correlations. 


(8) 


3 Results 


3.1 Spatiotemporal characteristics of water yield 


3.1.1 Temporal variation 

P, AET, PET, and simulated water yield in the WRB from 1985 to 2019 were shown in Figure 3. 
Annual average P over the study period was 542.09 mm with no significant temporal trend 
(Z=1.28, P>0.05). The minimum and maximum P of 364.26 and 766.12 mm occurred in 1997 and 
2003, respectively. The average PET was 695.55 mm, showing no significant trend (Z=0.43, 
P>0.05). The average AET was 458.81 mm, and showed a significant increasing trend (Z=1.93, 
P<0.05), suggesting that more than 80.00% of precipitation in the WRB was returned to the 
atmosphere in the form of evaporation every year. The water yield of the WRB fluctuated, and 
was correlated with precipitation. The maximum and minimum water yields of 261.62 and 6.08 
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Fig. 3 Annual precipitation, actual evapotranspiration (AET), potential evapotranspiration (PET) and water 
yield of the Weihe River Basin, China during 1985—2019 


202205.00019v1 


chinaXiv 


ChinaXivA (ERAT! 


WU Changxue et al.: Application of the InVEST model for assessing water yield... 


mm were found in 2003 and 1997, respectively. The annual average water yield was 83.14 mm, 
accounting for about 15.33% of annual precipitation. There was an insignificantly increasing 
trend in water yield over the entire study period (Z=0.94, P>0.05) at a rate of 4.84 mm/10 a, with 
a turning point occurring in 2003. 

3.1.2 Spatial variation 

As shown in Figure 4, the present study examined the spatial distributions of simulated annual 
water yield. For example, in 1985, the overall water yield of the basin was high. The main water 
producing areas were located in the Beiluo River Basin, the north Jinghe River Basin, and the 
upper reaches of the Weihe River. The annual water yield from 2005 to 2015 decreased from east 
to west and south to north. The main water-producing areas included the main stream of the 
Weihe River, and southern part of the Qinling Mountains. The low water yield areas located in the 
Jinghe River and Beiluo River. The global Moran's J obtained by the Geoda model ranged from 
0.38 to 0.68, indicating that water yield had a significant spatial agglomeration effect. 


(a) 1985 (b) 2005 (c) 2015 (d) Average 
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, m 
me <50 EE 50-100 EE 100-200 mm 200-300 mm >300 — Sub-basin boundary u 


Fig. 4 Spatial distribution of water yield in the Weihe River Basin, China. (a), 1985; (b), 2005; (c), 2015; (d), 
average. 


3.2 Response of water yield to P and LULC 


According to the principle of water yield model, P and AET were the main factors affecting the 
simulation results. However, AET was largely influenced by LULC. In this study, two scenarios 
were designed to explore the changes in water yield under different LULC and P conditions 
(Table 3). Scenario 1 was established to study the change in water yield under different P 
intensities, in which 1997, 2003, and 2010 represented dry, wet, and average years, respectively. 
Scenario 2 was established to study the change in water yield under different LULC conditions. 


Table 3 P and LULC under different scenarios 


Index Scenario 1 Scenario 2 
P 1997 2003 2007 1990 1990 1990 
LULC 2020 2020 2020 1990 2010 2020 


Note: P, precipitation; LULC, land use/land cover. 


3.2.1 Change in water yield under scenario 1 

The results of scenario 1 showed that water yield of the basin has changed significantly under 
different P conditions (Fig. 5). As P changed from less to more, the water yield of the basin 
changed accordingly. The maximum water yield during a wet year (675.40 mm) exceeded that 
during a dry year (133.53 mm) by a factor of 5, whereas it exceeded that during a normal year 
(365.20 mm) by a factor of 2. The average water yield of the basin varied greatly with changing P 
in 1997, 2003, and 2010 with 6.36, 264.63, and 72.96 mm, respectively. The spatial distribution 
of water yield was consistent with that of P. P in the study area was mostly concentrated in the 
southern and eastern areas, corresponding with the main water-producing areas of the WRB. 
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Fig. 5 Spatial distributions of annual precipitation (P, a-c) and water yield (d—f) in the Weihe River Basin, 
China under scenario 1 


3.2.2 Spatial correlation between water yield and P 


The bivariate Moran's J of P and water yield obtained by the Geoda model ranged from 0.30 to 
0.80, and showed a positive correlation (P>0.05). The LISA agglomeration maps of P and water 
yield showed aggregations of mainly high-high and low-low correlation relationships (Fig. 6). 
Areas with high-high relationships were mainly concentrated in the southern and eastern WRB. 
Areas with low-low relationships were mainly concentrated in northern Jinghe River and Beiluo 
River. 


(a) 1997 (b) 2003 (c) 2007 


Not significant W High-high EE Low-low E High-low Œ Low-high 0 100km 
— Sub-basin boundary 


Fig. 6 Spatial distribution of the relationship between annual precipitation and water yield in the Weihe River 
Basin, China in 1997 (a), 2003 (b) and 2007 (c) 


3.2.3 Change in water yield under scenario 2 


The results of scenario 2 showed that LULC had no significant effect on water yield, but had an 
impact on the spatial distribution of water yield (Fig. 7). There were little differences between the 
maximum and average water yield of the basin in 1990, 2010, and 2019 with 468.10, 466.88, and 
474.54 mm, respectively, and with 140.07, 140.9, and 142.85 mm, respectively. The distribution 
of land use types has a certain impact on the spatial characteristics of water yield. As shown in 
Figure 7, areas of higher water yield were correlated with grassland and farmland, such as 
northern Jinghe River and the upper reaches of the Weihe River. 


3.2.4 Water yield under different LULC 


The spatial analysis of ArcGIS10.2 was used to identify zones of water yield for each land use 
type under scenario 2. As shown in Figure 8, the fluctuations of water yield were different among 
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different LULC. The water yields of farmland, forestland, grassland, water body, and construction 
land continued to increase over the study period, whereas those of unused land first increased and 
then decreased. The maximum and minimum average water yield obtained for unused land and 
water body were 162.40 and 131.58 mm, respectively. 


C0 Farmland mm Forestland mi Grassland ™) Water body sm Construction land Unused land 


(a) 1990 (b) 2010 


(c) 2019 


(d) 1990 (e) 2010 () 2019 

Water yield (mm) Water yield (mm) Water yield (mm) 
468.10 466.88 

= 0.00 = 0.00 


0 100km 


Fig. 7 Spatial distributions of LULC (land use/land cover) types (a—c) and water yield (d—f) in the Weihe River 
Basin, China under scenario 2 


LULC of the WRB has changed significantly during the study period. As shown in Table 4, the 
areas of farmland, grassland, and construction land have changed significantly. Farmland area 
decreased by 4146.83 km?, whereas water yield increased by 20.00 mm. The areas of forestland, 
grassland and construction land increased by 498.42, 1624.8, and 1791.96 km?, respectively, 
whereas water yield increased by 11.34, 24.31, and 12.53 mm, respectively. The area of unused 
land decreased by 9.71 km?, whereas water yield decreased by 17.44 mm. The results indicated 
that the areas of forestland, grassland, unused land, and water body had positive correlations with 
water yield, whereas that of farmland had a negative correlation with water yield. However, 
grassland and farmland were the dominant land use types in the WRB, accounting for more than 
80% of the basin's area in 2020. Therefore, water yield was larger in grassland and agricultural 
land. 


Table 4 Changes of LULC (land use/land cover) types of the Weihe River Basin China from 1990 to 2020 


Construction Unused ‘Water 


fides pri Poy re land land body ne 
km’) (km’) (km’) 

Farmland 32,855.50 3589.30 18,335.17 3383.78 68.70 491.49 58,723.94 
Forestland 2933.04 11,916.83 6314.91 177.19 30.31 60.66 21,432.95 
Grassland 16,536.86 6219.80 25,866.74 703.06 58.78 207.40 49,592.64 
Construction land 1770.24 117.94 448.83 741.28 4.94 43.82 3127.04 
Unused land 61.76 39.77 60.20 7.99 36.59 3.86 210.16 
Water body 419.71 47.73 191.59 85.70 1.14 128.26 874.14 


Total 54,577.11 21,931.37 51,217.44 5099.00 200.45 935.50 133,960.87 


202205.00019v1 


chinaXiv 


ChinaXiv@ ERAT! 


JOURNAL OF ARID LAND 
200 5 
me 1990 3 2010 wm 2019 == Average 
150 4 
100 4 


Water yield (mm) 


50 4 


Farmland Forestland Grassland Water body Construction land Unused land 
Land use type 


Fig. 8 Water yield of different LULC (land use/land cover) types in the Weihe River Basin, China 


4 Discussion 


4.1 Effects of various factors on water yield 


A combined examination of the results of trend and spatial analysis as well as scenario analysis 
showed that precipitation was the most direct factor affecting regional water yield, which is 
consistent with the relevant research results (Jiang et al., 2016; Li et al., 2021). As for AET, it 
directly participates in the hydrological cycle process (Lewis and Allen, 2017). With emergence 
of challenges resulting from global climate change, AET has gradually increased in recent 
decades (Fig. 3). The increase in AET may accelerate the hydrological cycle process, and affect 
the temporal and spatial distributions of hydrological elements (Gusev et al., 2019; Cheng and Li, 
2020). AET is not only affected by meteorological factors (temperature, wind speed, relative 
humidity and sunshine hours), but also by the land use conditions (Lang et al., 2017). 

There is variation in water yield among different land use types due to differences in soil water 
content, evapotranspiration capacity, litter water holding capacity, and canopy interception. 
Unused land had the largest water yield because a greater proportion of precipitation directly 
penetrated into the ground or formed surface runoff (Lang et al., 2017). The increase in 
construction land in the basin has not only increased the demand for water resources, but also 
changed the underlying surface conditions. Construction practices transform the ground surface 
into an impermeable layer by removing the vegetation. These conditions lead to a decrease in the 
evapotranspiration and a higher water yield capacity (Sterling et al., 2013; Anache et al., 2017). 
The water yield capacity of water body is weak due to strong evaporation on the water surface. 
The water yield generated by forestland was relatively lower because of high transpiration and 
interception of water by the deep root system, litter layer, and dense canopy (Vose et al., 2011; Li 
et al., 2021). The effect of grassland and farmland on the redistribution of precipitation was 
similar to that of forestland. However, the regulation effect of grassland and farmland was 
relatively weak due to a lower canopy coverage and shallower root depth. In addition, the wide 
distributions of these two types of lands in the basin has resulted in their dominant contributions 
to water yield. Previous studies have shown that grassland was the optimal land use patterns for 
maintaining hydrology (Li et al., 2021). Although there have been considerable changes to the 
land use in the basin during the study period, these changes had no significant effect on water 
yield. This result can be attributed to the mutual transformation of various land types offsetting 
any dominant trend in water yield change due to land use changes (Nie et al., 2011). 

Land use is the product of human activities, emphasizing human use of natural lands (Lambin 
et al., 2003; Goldewijk and Ramankutty, 2004; Liu et al., 2013), which can affect processes and 
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components of ecosystems. A series of soil and water conservation projects have been conducted 
in the basin since the 1950s (Mu et al., 2007). These include the construction of terraces and 
sediment dams, afforestation, plant restoration, pasture improvement, and the returning of 
farmland to forestland or grassland. These projects have significantly improved the ecological 
environment and water supply functions of the basin (Wang et al., 2017; He et al., 2021). The 
principle of water balance indicates that the difference between water yield and measured runoff 
is the water consumed by agricultural, urban and industrial activities, and the changes in reservoir 
storage (Li et al., 2020). The current study calculated the runoff of the WRB as sum of measured 
runoff of the Huaxian and Zhuangtou hydrological stations (Fig. 1) over years. The results 
showed an upward trend in annual water consumed (Z=1.56, P<0.01), which has increased 
pressure on water resources. 


4.2 Limitations of the nVEST model 


The InVEST model has been broadly used to evaluate ecosystem service functions, and has 
achieved good results (Lang et al., 2017; Kim and Jung, 2020; Daneshi et al., 2021). However, 
there are some uncertainties in the model related to model setup and its simplified algorithm 
(Sharp et al., 2020). For example, the model does not consider complex terrain factors, and is not 
able to rigorously describe the water balance process under complex underlying surface 
conditions (Jiang et al., 2016). Topography affects the climate by changing regional hydrothermal 
conditions. These changes influence the growth and structure of vegetation, the accumulation of 
litter and the physicochemical properties of soil, which consequently affect water yield (Jia et al., 
2014; Maurya et al., 2016). 

In addition, the present study obtained root depth and soil data for input into the model from 
the global soil database. The low spatial resolution of these data affected the simulation accuracy 
of the model to a certain extent. Moreover, the value of seasonal parameter Z used in the present 
study was slightly different from that applied in similar watersheds. This indicates that differences 
in natural conditions will lead to large changes in Z, even between similar basins. Therefore, the 
value of Z is important to verify before the INVEST model is applied. 


5 Conclusions 


The Weihe River is the largest tributary of the Yellow River, and has strategic significance for the 
protection of ecological environment and water resources management in Northwest China. The 
present study applied the nVEST model to quantitatively evaluate the water yield of the WRB 
from 1985 to 2019. In addition, the response of water yield to climate factors and land use types 
was explored. 

The average annual water yield of the study area was 83.14 mm, and yearly water yield over 
the study period showed a slight increasing trend. The main stem of the Weihe River and southern 
region were the main water-producing areas. Water yield was comprehensively affected by 
climate and land use factors. Precipitation was the most direct influencing factor, with the spatial 
distribution of precipitation corresponding with that of water yield. In contrast, land use did not 
have a significant effect on water yield. However, there was variation in water yield among 
different land use types. Unused land had the largest water yield capacity, but farmland and 
grassland contributed the most to the total water yield of the basin due to their wide distribution. 
The results of the present study can act as a reference for formulating reasonable and efficient 
water resources allocation schemes. In addition, the successful application of the InVEST model 
in the current study can guide its application in related researches under similar natural 
conditions. 
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